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The total electron content Nt and the equivalent slab thickness r of the midlatitude ionosphere 
are determined from differential Doppler measurements on the radio transmissions from the Transit 
4A satellite and the ionosonde data near Washington, D.C., and Ottawa. The trends of diurnal, 
seasonal, and sunspot cycle variations of Nt and r at midday are described. It is found that in 1962 
both Nt and r were considerably decreased in magnitude as compared to the corresponding values near 
the peak of the sunspot cycle. The inferred mean electron-ion temperature is found to be generally 
in excess of the neutral gas temperature. The ratio Te/T,, in the daytime is estimated to be 1.8 in 
summer and about 1.1 to 1.2 in winter, which may reflect corresponding changes in the ratio Q/N^, 
where Q is the heat input to the electrons and Ne is the electron density. The relationship between 
T and Kp appears to be rather ambiguous in contrast with the positive correlation between the neutral 
gas temperature and Kp. 



1. Introduction 

The radio transmissions from artificial earth satel- 
lites have been very useful in the exploration of the 
ionosphere, particularly for the measurements of the 
total electron content, Nt, up to the height of the satel- 
lite and for studying ionospheric irregularities. The 
characteristics of the satellite radio signals such as 
amplitude, frequency, phase, and direction of arrival 
may all undergo significant changes during their pass- 
age through the ionosphere, provided the frequency 
of the radiation is not too large in comparison with 
the critical frequency of the F region of the ionosphere. 

Measurements of the Faraday fading of the satel- 
lite signals have been used to calculate Nt by [Garriott, 
1960; Little and Lawrence, 1960; Blackband, 1960; 
Yeh and Swenson, 1961; Blumle, 1962; Lawrence et al., 
1963; Rogers, 1964] and others. Measurements 
involving Doppler frequency shifts have been used to 
calculate Nt by [Aitchison et al., 1959; Ross, 1960; 
de Mendonca, 1962, and Hibberd 1964]. Garriott 
and de Mendonca [1963] showed that the most accurate 
method of obtaining Nt is a hybrid analysis using simul- 
taneous Faraday and Doppler data. Recently, Tithe- 
ridge [1964] determined the diurnal and seasonal 
variations of Nt and the scale height, H, of the F 
region, by measuring the vertical angle of arrival of 
the 20-MHz radio signals from Explorer 7 at Auckland 
over a period of 16 months in 1960-61 (refraction 
method). Measurements of the equivalent slab 
thickness, r, defined as the ratio of Nt measured by 
satellite observation to the corresponding maximum 
electron density of the F region, A^max? have been 
discussed by [Ross, 1960; Ross and Anderson, 1962; 
Rogers, 1964; and Hibberd, 1964]. 



The present paper describes the results of the 
analysis of the differential Doppler data received near 
Washington, D.C., by the Johns Hopkins University 
on two pairs of harmonically related frequencies (54- 
324 MHz and 150-400 MHz) transmitted from Transit 
4A during 1962, a period of low solar activity. The 
method of calculating Nt has been described elsewhere 
[sec. 3.4, Garriott and de Mendonca, 1963]. Mean 
diurnal, seasonal, and solar cycle variations of Nt and 
T of the midlatitude ionosphere will be described and 
compared with the results obtained by other workers. 

2. Experimental Data 

Of 702 sateUite recordings available from February 
1962 to September 1962, only 259 were used in the 
present study. The remaining 443 recordings were 
rejected either because the satellite ephemeris was 
not available or the recordings were so much disturbed 
that the "proximal point" of the satellite could not be 
determined. The intersection of the propagation ray 
path with a concentric spherical shell at the height of 
350 km is called the ionospheric point, and its vertical 
projection on the ground is the subionospheric point. 
A few sample plots of A^r against geographic latitude of 
the subionospheric point are shown in figure 1. While 
it is easy to obtain an almost continuous distribution of 
Nt over a small range of latitudes about the satellite 
receiving station, it is difficult to obtain the correspond- 
ing distribution of r owing to the sparsity of ionosonde 
stations. Therefore, we shall restrict ourselves to the 
diurnal and seasonal variations of the ionosphere over 
two midlatitude stations, namely, Washington, D.C. 
(40° geographic) and Ottawa (45°.4 geographic) for 
which both Nt and A^max values were available. A^max 
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Figure 1. Some sample plots o/ Nt against geographic latitude of the subionospheric point 

of Transit 4A satellite in 1962. 

To represents the proximal point of satellite passage. 

930 



values for Washington and Ottawa were taken from the 
published ionosonde data. All A^max values were inter- 
polated to the corresponding local time of the satel- 
lite passage. 

3. Results of the Analysis 
3.1. Diurnal and Seasonal Variations of iVr 

Figures 2a and 2b show the values of Nt up to the 
height of the satelhte (usually greater than 880 km) 
over Washington and Ottawa respectively, obtained 
from the analysis of the differential Doppler data of 
Transit 4A satellite during the period February 1962 
to September 1962. These values of Nt were scaled 
from the latitudinal plots of Nt for both the northbound 
and southbound passages of the satellite. The local 
times and the months of observations are shown at 



the bottom of the figures. The smoothed sunspot 
number, R, decreased from 42 to 32 during the ob- 
serving period. Although 13 geomagnetic disturb- 
ances were reported during this period, none of them 
was classified as severe. The storm indices ranged 
from ^ = 10 to ^ = 40, with a median of about 25. 
(Storms with A indices of about 50 or higher are classi- 
fied as severe.) Therefore, all the available values 
of Nt were plotted in figures 2a and 2b. Although 
there is scatter in the values of Ar on the order of ± 30 
percent the diurnal variation of about 6 to 1 can be 
clearly seen both at Washington and Ottawa. Scatter 
of Nt is believed to occur mainly due to large day-to-day 
changes in the ionosphere and is not an indication of 
measurement error. 

For the purpose of comparing Nt near the latitude 
of Washington (40°N) and that of Ottawa (45°N), the 
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Figure 2. The total electron content Nj of the ionosphere deduced from Transit 4A differential Doppler 
measurements at (a) 40°N and (b) 45°N geographic latitudes. 
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Figure 3. Mean diurnal variation of Nt at 40°N and 45°N (geographic). 



mean diurnal curves drawn manually through the 
pomts in figures 2a and 2b are shown in figure 3. The 
diurnal variations of Nt at these stations appear to 
be similar with a daytime peak that occurs in the 
afternoon and a minimum before sunrise. Noon 
values of Nt at Washington, however, were somewhat 
higher than those at Ottawa. This is not unexpected 
since a similar trend is found in the noon A^max values 
at these stations. It should be noted that the noon 
values of Nt in March-April 1962 were slightly higher 
than the corresponding values in June-July 1962. 
This seems to be due to seasonal variation of Nt-, 
although a small part of this seasonal variation may 
be attributed to the declining solar activity. During 
1962, both the diurnal and seasonal variations of Nt 
at midlatitudes diminished in amplitude as compared 
to the variations observed by Lawrence et al., [1963] 
for the Boulder and Stanford data obtained from the 
Faraday rotation of the signals from Sputnik 3 in 
1958-59. They observed a seasonal variation in the 
noon values of Nt of about 2 to 1 from winter to summer 
and a diurnal variation of about 10 to 1. 



3.2. Mean Diurnal and Seasonal Variation of the 
Equivalent Slab Thickness, r 

The equivalent slab thickness of the ionosphere, 

Nt 
^^N~~ 

is evaluated from the satellite observations coupled 
with the corresponding values of A^max from ionosonde 
observations [Ross, I960]. For a Chapman-alpha 
electron density profile, we may deduce that A^t = 4.13 
H A^max? or that t = 4.13 f/ where H is the scale height 
of the ionizable constituent. Thus, the measurement 



of T can be easily referred to the scale height (or tem- 
perature) of the ionized layer, if the electron density 
profile is reasonably close to a Chapman model. At 
night, this has been shown to be expected [Duncan, 
1956; Dungey, 1956] and, for the daytime, Bowhill 
[1962] and Nisbet [1963] have obtained equilibrium 
solutions which have the same electron content as 
the Chapman layer, to within 10 percent. In the 
absence of thermal equilibrium, the implication of the 
slab thickness is not so clear, but it should represent 
at least approximately the mean electron-ion tempera- 
ture of the plasma. 

Figure 4 shows the mean diurnal variation of r 
over Ottawa for the period February to September 
1962 calculated in two different ways. The curve 
represented by (X) has been determined from the 
smoothed values of A^r and A^max at 1-hr intervals from 
to 23 hr local mean time. The second curve of r 
shown by (0) has been determined by grouping the 
individually calculated values in similar 1-hr intervals. 
There appears to be good qualitative agreement 
between the two diurnal curves of r. It is clear that r 
exhibits a significant diurnal variation at least of the 
order of 30 percent, with a peak near sunrise period. 
The actual diurnal variation may not be quite as large 
as seen in figure 4, since this mean diurnal curve has 
been developed from the data extending over a period 
of eight months. The occurrence of a peak in r near 
sunrise has been reported by many investigators [Evans 
and Taylor, 1961; Lawrence et al., 1963; Titheridge, 
1964]. This has been explained by the fact that 
sunrise occurs earlier at greater heights so that the 
ionization above the peak of F layer will have con- 
siderably increased by the time the electron density 
at the peak begins to increase [Titheridge, 1964]. 
Also, measurements of electron temperature Te, 
made by using incoherent backscatter technique 
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Figure 4. Mean diurnal variation of the equivalent slab thickness 
7 over Ottawa in 1962, calculated from the mean diurnal variation 
of Nt and N^ax W o,f^d from individual observation of Nt and 

Nn,ax (0). 



[Bowles et al., 1962] and by the Ariel satellite [Will- 
more et al., 1962] show evidence of a significant 
departure from thermal equilibrium between the elec- 
trons and ions near sunrise. This would increase the 
mean temperature of the electron-ion gas resulting in 
greater slab thickness at sunrise. 

Assuming a Chapman profile for the daytime F 
region and atomic oxygen of mass unit =16 AMU, 
the noon values of r at Ottawa in 1962 yield a scale 
height of about 87 km which corresponds to a tempera- 
ture of about 1500 °K. This is in excess of the neutral 
gas temperature Tg, which is about 1000 °K as calcu- 
lated by Nicolet [1963] from satellite drag analyses. 
This temperature difference may be attributed to the 
absence of the thermal equilibrium between charged 
and neutral particles. The change in r from day to the 
early part of the night corresponds to a temperature 
change of about 470 °K. It is found that in 1962, 
the inferred mean electron-ion temperature in the early 
part of the night is greater than the corresponding 
neutral gas temperature by about 150 °K. Thus, there 
is a definite tendency towards thermal equilibrium 
between the charged and neutral particles in the night- 
time ionosphere. 

Seasonal variation of the noon values of r for mid- 
latitude ionosphere have been measured by different 
workers. Their results, converted to r where neces- 
sary, are summarized in table 1. 

An examination of table 1 indicates that the seasonal 
variation of r at noon in the northern hemisphere 
varies between 30 to 50 percent with a mean of about 
40 percent, while in the southern hemisphere it is only 
about 20 percent. One possible reason for such wide 
variations in the noon values of r may be the latitudinal 
differences in A^max at different stations, particularly 
in view of the fact that the mean electron-ion tempera- 
ture depends critically on the ratio Q/Nj, where Q 
is the heat input to the electrons due to solar ionizing 
radiation and A^e is the electron density [Hanson, 
1963]. 



Table 1. Seasonal variation ofr at noon 



Authors 


Year 


Place 


Geomag. 
lat. 


Geomag. 
long. 


Ratio of T 
in summer 
and winter 

(Ts/Tw) 


Lawrence 
et al. 


1958-59 


Stanford 
and 
Boulder 


43. yN 
48.9°N 


298.4° 
316.4° 


1.30 


Ross and 

Anderson 


1959 


University 
Park, Pa. 


.52.2''N 


349.1° 


1.51 


Rogers 


1959 


Jodrell 
Bank 


.56.:rN 


82.5° 


1.49 


Hibberd 


1960-61 


University 
Park, Pa. 


.52.2°N 


349.1° 


1.39 


Munro 


1959 


Baxland 


42.4°S 


226.5° 


1.16 


Titheridge 


1960-61 


Auckland 


41.4°S 


252.7° 


1.20 



3.3. Equivalent Slab Thickness (r) of the Ionosphere 

Over Washington and Ottawa 

The mean diurnal curves of r for Washington (40°N) 
and Ottawa (45°N) are shown in figure 5 for comparison. 
It is observed that r at Ottawa is systematically higher 
than T at Washington for much of the time except 
around midnight hours. At midday, the difference in 
T at these stations is about 30 km. This would imply 
that the mean temperature of the electron-ion gas in 
the ionosphere over Ottawa must have been about 
125 °K greater than that over Washington. This con- 
clusion seems to be consistent with the electron tem- 
peratures observed with Aerial satellite which showed 
that Te increased with increasing latitude. Similar 
latitudinal trend in [{Te + 7^0/2] at 500 km was noticed 
in the analysis of the Alouette data in the fall of 1962 
[Bauer and Blumle, 1964]. 

3.4. Solar Cycle Variation of the Midday Total 

Electron Content {Nt) at Midlatitudes 

It is a well-known fact that the critical frequency, 
/oF2, of the F2 layer varies markedly with solar cycle. 
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Figure 5. Comparison of mean diurnal variation of r over Wash- 
ington, D.C, and Ottawa in 1962. 
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Figure 6. Solar cycle variation of Nj of midlatitude ionosphere 
and the mean sunspot number, R. 



being greatest at sunspot maximum. The worldwide 
results have been summarized by the expression 
(/o/^2)2oc(l + 0.02 R), where R, is the mean Zurich 
sunspot number [Ratcliffe and Weekes, I960]. It is 
of interest to determine a corresponding relationship 
between Nt and R. In figure 6 are plotted the mean 
midday values of Nt at midlatitudes against R in 
different months measured by various workers during 
the period of declining solar activity. It is evident 
that Nr increases linearly with /?, when the data are 
grouped by seasons, and we may express this empirical 
relationship as 

NT = a[l-^- (i?-40)] X W elec/m2 

for R ^ 40. The parameters a and b are constants 
which are different for different seasons. The vari- 
ation of a, b, and the ratio b/a are shown in table 2. 



Table 2. Seasonal variation of 

a, b, and - for R^40 at 

a -^ 

midlatitudes 



Season 


a 


b 


bla 


Summer 


1.5 


0.017 


0.011 


Winter 


1.6 


.036 


.022 


Equinox 


2.0 


.05 


.025 



The ratio bj_a represents the 'sensitivity' of Nr for 
changes in R. Examination of table 2 shows that the 
ratio bla seems to be comparable in magnitude in 
winter and equinox but appreciably smaller in summer. 
Thus, it appears that, in winter and equinox, Nt 
responds much the same way to the changes in solar 
activity, even though the absolute values of Nt might 
differ in these seasons. However, it is. clear that in 
summer, Nt increases with solar activity less rapidly 
than either in equinox or winter. This explains the 
enhanced seasonal variation of midday Nr during a 
period of high solar activity. Taylor [1963] and more 
recently Hibberd [1964] combined their mean midday 
values of Nt for summer and winter seasons with those 
of others and showed that, in winter, Nt increased 
linearly with 10.7 cm solar radio flux more rapidly 
than in summer, which is in accord with the present 
results. 

3.5. Solar Cycle Variation of the Midday Slab 
Thickness r at Midlatitudes 

It is generally accepted that solar extreme ultra- 
violet radiation is an important source responsible for 
heating of the F region [Hunt and van Zandt, 1961]. 
Satellite drag measurements show that the thermo- 
pause temperature and hence the neutral particle 
densities at F region heights vary linearly with the solar 
activity [Nicolet, 1963]. In figure 7 are plotted the aver- 
age midday values of r_for each season against the 
mean sunspot number R^ The empirical linear rela- 
tionship between r and R may be written as 

T (km) = 270 (1 -h 0.005 R) for summer, 

T (km) = 240 (1 -h 0.005 R) for equinox, 

and T (km) = 210 (1 + 0.005 R) for winter. 

Clearly, r increases progressively from winter to 
summer. The seasonal variation of r has been 
observed to be about 30 percent during the period 
of declining solar activity from 1958 to 1962. Addi- 
tional scales for scale height H and mean electron-ion 
temperature [{Te + Ti)l2)], assuming mean ionic mass 
equal to 16 AMU, are also shown in the figure 7. 
Since Nt ^ A^max ^ ^^ ^^ percent of the solar cycle 
variation of Nt, in equinox and winter, can be attri- 
buted to A^max variation and 25 percent to the tempera- 
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Figure 7. Solar cycle variation of the equivalent slab thickness 
of the midlatitude ionosphere and the mean sunspot number, R. 



ture variation, whereas, in summer, A^max and 
temperature seem to make approximately equal con- 
tributions to Nt variations. 

Ross and Anderson [1962] and Hibberd [1964] re- 
ported positive correlation between daytime r and the 
solar activity as measured by either decimetric 
solar radio flux or the mean sunspot number. For 
southern winter, Titheridge [1964] determined the 
variation of H with the mean sunspot number given by 

//= 54(1 + 0.0023^) km, 

which is approximately half as large as the correspond- 
ing variation in the northern hemisphere. 

The mean electron-ion temperature inferred from 
the measurements of r seems to be higher than the 
corresponding thermopause temperature calculated 
by Nicolet [1963] in all seasons. The mean electron- 
ion temperature (Te-\-Ti)/2 at midday inferred from 
Sputnik 3 data and Transit 4A data, assuming a 
Chapman profile for the F region, are shown in table 
3, along with the neutral gas temperature (Ty) deter- 
mined by Nicolet. Assuming Tj — Tg, the ratio of 
TelTff has been calculated for different seasons. It is 
seen that Te/Tg is of the order of 1.8 in summer but 
tends to be progressively lower in equinox and winter. 
This seasonal behavior of the ratio Te/Tg reflects 
similar variations in the ratio of Q/Nj where Q is 
the heat input to the electrons and A^e is the electron 
density. 



Figure 3. Seasonal variation of Tp and Tg at 
Midday 

















Year 


Season 


K 


n-Ti 


{n + Ti)l2 


Te 


TelT, 


1958-59 


Wintt-r 


180 


1670 


1720 


1770 


1.06 




Etjiiinox 


170 


1580 


1800 


2020 


1.28 




Summer 


155 


1500 


2100 


2700 


1.8 


1962 


Winter 


50 


1000 


1125 


1250 


1.25 




Equinox 


40 


975 


1220 


1465 


1.5 




Summer 


35 


950 


1337 


1724 


1.82 



3.6. Scale Height and Magnetic Activity 

As mentioned earlier, no severe geomagnetic dis- 
turbances occurred during the observing period in 
1962 for which the data were available. However, 
Hibberd [1964] showed that, at least during some in- 
tense disturbances, the scale height increased sub- 
stantially, in agreement with the results of Lawrence 
et al., [1963] but not with Titheridge [1964] who 
observed that the scale height did not vary by more 
than a few percent with magnetic activity. In order 
to further investigate this question, we reexamined the 
Sputnik 3 data of Lawrence et al., [1963] for possible 
storm effects. A separate 'quiet' day curve of the 
scale height for each season was determined from the 
mean diurnal variation curve obtained by these 
authors, after correcting it for seasonal and solar 
cycle variation. An assumption was made that the 
form of the diurnal variation was independent of the 
season. The scale height values for Kp^ 4' were 
regarded as 'disturbed' and these were found to be 
higher as well as lower than the corresponding 'quiet' 
day values. The effect of Kp on H is summarized in 
table 4. 



Table 4. Effect of Kp on H for Kp ^ 4 in 1959 



Season 


Nu 


mber of cases 


show 


ng 




Increase in // 


Decrease 


in// 


Summer 




9 




17 




Equinox 




15 




6 




Winter 




6 




5 





It is suggested by table 4 that increases in H with 
magnetic activity are more common in equinox while 
the opposite seems to be true in summer. In winter, 
there is no significant tendency for H to increase with 
magnetic activity. Thus, at a midlatitude station, the 
correlation between the scale height or mean electron- 
ion temperature appears to be rather ambiguous in 
contrast with the definite positive correlation between 
the neutral gas temperature and the magnetic activity 
[Jacchia, 1964]. The reason for this apparent dis- 
crepancy may be related to the nature of different heat- 
ing mechanisms for the neutral and electron-ion gas. 
The heating of the neutral atmosphere may be accom- 
plished through processes such as hydromagnetic wave 
absorption or joule losses which do not directly con- 
tribute to the plasma temperature. However, the solar 
EUV radiation is generally assumed to be the main 
source of heat in the quiet midlatitude daytime iono- 
sphere. Ionospheric storms are considered to be 
either 'positive' or 'negative' depending upon whether 
^max increases or decreases during the disturbance 
[Maeda and Sato, 1959]. For a given value of heat 
input Q, which may be relatively constant in a storm, 
these variations of A^max niay alter the ratio QIN^, 
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which determines the temperature difference 
(Te — Ti), in such a way as to cause an increase or 
decrease in the mean electron-ion temperature and 
therefore the scale height. 

4. Conclusion 

The total electron content, Nt, and the equivalent 
slab thickness, r, of the midlatitude ionosphere has 
been determined from the differential Doppler data 
of Transit 4A for the period February to September 
1962. Significant diurnal and seasonal variations in 
Nt and r were observed. It was found that the ampli- 
tudes of these variations were diminished considerably 
as compared to the earlier observations made with 
Sputnik 3 near the peak of the sunspot cycle. The 
mean diurnal curve of r showed a peak near sunrise. 
The change in r from day to night was of the order of 
30 percent which implied a fall in temperature of the 
electron-ion gas of about 470 °K. The seasonal vari- 
ation of midday values of r from summer to winter in 
the northern hemisphere ranged from about 30 to 50 
percent with a mean of about 40 percent. This is 
about twice the seasonal variation of r observed in 
the southern hemisphere. The effect of solar activity 
on the daytime equivalent slab thickness, r could be 
summarized by the expression of the form r oc (l-h0.005 
R) where R is the mean sunspot number. The daytime 
T over Ottawa was about 30 km greater than that over 
Washington which implied that the mean electron-ion 
temperature in the F region over Ottawa was greater 
by about 125 °K than that over Washington. 

For R ^ 40, it was found that the solar cycle depend- 
ence of noon values of Nt could be expressed as 

A^r = [a + 6 (^ - 40) ] X W elec/m2 

where a and b are constants which are different for 
different seasons. 

The mean electron-ion temperature in the daytime 
ionosphere was found to be generally in excess of the 
corresponding neutral gas temperature determined 
from the satellite drag measurements. This is inter- 
preted as an indication of the absence of the thermal 
equilibrium between charged and neutral particles. 
This inference is based, however, on an assumed pro- 
file shape for the electron distribution. Assuming 
Ti — Tfj, the daytime ratio of Te/Tg was estimated to be 
about 1.8 in summer and about 1.1 to 1.2 in winter. 
Reexamination of Sputnik 3 data for 1958-59 showed 
that the relationship between the scale height or the 
mean electron-ion temperature and magnetic activity 
appears to be rather ambiguous in contrast with the 
definite positive correlation between the neutral gas 
temperature and magnetic activity. 



The recordings of Transit 4A signals in 1962 have 
been very generously loaned to Stanford University 
through the courtesy of Guy Worsley and George 
Weiffenbach from the Applied Physics Laboratory of 
the Johns Hopkins University under Contract No. 



w62 -0604c with the Department of the Navy. The 
data were originally recorded in connection with 
the tracking of Transit sateUites under the direction 
of APL. Fernando de Mendonca arranged for this 
loan while he was at Stanford University, and much 
of the data reduction has been supervised by S. C. 
HaU. The analysis has been supported by NASA 
Grant NsG 30—60. R. V. Bhonsle is on leave of 
absence from the Physical Research Laboratory, 
Ahmedabad, India. 
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